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SUMMARY 


Theoretical aspects of a new capability to determine the 
vibratory response of turbosystems subjected to aerodynamic 
excitation are presented. 

Turbosystems such as advanced turbopropellers with highly 
swept blades, and axial-flow compressors and turbines can be 
analyzed using this capability. The capability has been 
developed and implemented in the April 1984 release of the 
general purpose finite element program NASTRAN. 

The dynamic response problem is addressed in terms of the 
normal modal coordinates of these tuned rotating cyclic 
structures. Both rigid and flexible hubs/disks are considered. 
Coriolis and centripetal accelerations, as well as differential 
stiffness effects are included. 

Generally non-uniform steady inflow fields and uniform flow 
fields arbitrarily inclined at small angles with respect to the 
axis of rotation of the turbosystem are considered as the sources 
of aerodynamic excitation. The spatial non-uniformities are 
considered to be small deviations from a principally uniform 
inflow. Subsonic and supersonic relative inflows are addressed, 
with provision for linearly interpolating transonic airloads . 

A stand-alone pre-processor program, independent of NASTRAN, 
has been additionally developed to generate the vibratory 
®i^losds on the blades of these turbosystems . 



Both NASTRAN and pre-processor capabilities are operational on 
the CRAY 1— S computer system at NASA ' s Lewis Research Center . 


The work was conducted under Contract NAS3-24387 from NASA 
LeRC, Cleveland, Ohio, with Mrs. Marsha Nall as the Program 
Monitor . 
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1. INTRODUCTION 


In a series of related efforts, NASA's Lewis Research Center 
has sponsored the development of a number of analytical 
capabilities addressing the static and dynamic problems of 
axial-flow turbomachines. To benefit from the state-of-the-art 
structural modelling and analyses techniques, these analytical 
developments have been implemented in the general purpose finite 
element program NASTRAN, These additional capabilities in 
NASTRAN, based on a unified approach to representing and 
integrating the structural and aerodynamic aspects of the 
tu rbomachinery problems , are listed below: 

1. Static aerothermoelastic design/analysis of axial-flow 
compressors (Refs . 1-3) , 

2. Modal flutter analysis of axial-flow turbomachines (Refs . 
1-3) , 

3. Forced vibration analysis of rotating cyclic structures , 
direct approach (Refs . 4,5) , and 

4. Modal flutter analysis of advanced turboprops with highly 
swept blades (Refs . 6,7). 

The new capability discussed in this report for modal forced 
vibration analysis of aerodynamically excited turbosystems , 

1 . extends the direct (as compared to modal) approach forced 
vibration analysis capability listed above to treat the 
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dynamic response problem in terms of natural modal 
coordinates of the rotating cyclic structure of the 
turbosystem, 

2. draws upon, and extends, the turboprop flutter analysis 
capability listed above for the generation of reactionary 
aerodynamic forces due to the vibratory motion of the 
turbosystem blades, and 

3. requires the specification of imposed oscillatory airloads on 
the turbosystem. (As per NASA's requirements, the generation 
of oscillatory airloads for the purposes of the present 
contract has been accomplished by the development of an 
independent pre-processor program discussed in detail in Ref. 
8 and briefly in Section 3 of this report . ) 

Mathematical description of the problem and its solution are 
presented in Section 2. 

Application examples on two advanced turboprops at various 
operating conditions are given in Section 4. 

Section 5 presents and discusses the results of Section 4 
examples . Conclusions and recommendations are reported in 
Sections 6 and 7 , respectively . 
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2 . MATHEMATICAL FORMULATIONS 


2.1 PROBLEM DESCRIPTION 

Figure 2.1 shows a single-rotation advanced turboprop, as an 
example of turbosystems, operating in a generally non-uniform 
steady inflow field. 

Although the absolute inflow field does not change with time, 
the rotation of the turboprop results in velocities with 
oscillatory components relative to the blades. Relative 
velocities with harmonic components at the rotational frequency 
also exist in uniform flow fields when the turboprop axis of 
rotation is misaligned with the absolute flow direction. 

Given such operating conditions , it is desired to, 

1 . determine the oscillatory loading distributions over the 
blades of the turboprop at various excitation frequencies , 
and 

2 . determine the resulting vibratory response (displacements , 
stresses , etc . ) of the turboprop. 


2.2 APPROACH 


Due to its eventual implementation as part of the Bladed 
Shrouded Disks Computer Program in NASTRAN at NASA LeRC, the 
approach adopted for the formulation and solution of the above 
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Spatially Non-uniform 



Figure 2 . 1 Advanced Turboprop in a 

Generally Non-Uniform Steady 
Inflow Field 


problem is derived from those described in Refs. 4 and 6. 

The following structural and aerodynamic aspects are 
considered in the theoretical development: 

Structural Aspects 

1. Structural representation 

- Structures of turbosystems are treated as tuned cyclic 
structures with identical masSf stiffness# damping# and 
constraint properties for all cyclic sectors. The 
structural modelling capabilities of NASTRAN for 
rotationally cyclic structures are admitted. 

2. Hub flexibility in relation to that of blades 

- A relatively rigid hub effectively decouples the 
structural behavior of the blades # any one of which can 
then be modelled as an independent structural member 
coupled only aerodynamically with the other blades of 
the turboprop. 

A hub of stiffness comparable to that of blades 
structurally couples the blades # in addition to their 
aerodynamic coupling . The structural behavior of the 
multi-bladed turboprop as a whole # as opposed to that of 
an individual blade # then influences the dynamic 
response of the turboprop. From a structural modelling 
viewpoint# however# with the use of analysis invoking 
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cyclic syimnetry features (Ref. 4), it still suffices to 
model only a one-bladed segment of the turboprop. 

3. Differential (or incremental) stiffness 

- The steady state loads, particularly centrifugal loads, 
acting on the relatively thin blades of advanced 
turbopropellers noticeably affect the geometry and 
effective stiffness of the blades. This, in turn, 
alters the relative inflow conditions at the blade 
leading edge, and also influences the natural 
frequencies of the turbopropeller blades at any 
operating condition. In general, the significance of 
the differential stiffness effects decreases with 
decreased rotational speed and increased elastic 
stiffness of the structure . 

4 . Coriolis and centripetal acceleration effects 

- These effects are present during vibrations of rotating 
turboprops , and are accounted for in the forced 
vibration analysis capability of the Bladed Disks 
Computer Program. 

5 . Use of normal modes of the turboprop in formulating and 
solving the dynamic response problem 

- The normal modes of tuned cyclic turbosystems can be 
grouped in mutually independent families corresponding 
to permissible circumferential harmonic indices . Hence , 
the total response of the turbosystem can be obtained by 
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combining the mutually independent responses for each of 
the applicable circumferential harmonics. 

Aerodynamic Aspects 

1. Aerodynamic modelling 

This is essentially dictated by the unsteady aerodynamic 
theories used to determine the unsteady blade loading 
distribution. The Bladed Disks Computer Program 
presently utilizes two-dimensional cascade aerodynamics 
in a strip theory manner. The blade model, accordingly, 
is a series of chordwise strips stacked spanwise to 
cover the entire blade surface as shown in (Figure 2.2) . 

2. Two-dimensional cascade aerodynamics for oscillatory blade 
motion 

Subsonic and supersonic aerodynamic routines , currently 
part of the Bladed Disks Computer Program, with 
modifications to account for blade sweep , if present , 
are used to generate the reactionary aerodynamic blade 
loading distributions . 

3. Two-dimens ional cascade aerodynamics for oscillatory relative 
inflow 

Versions of subsonic and supersonic routines capable of 
handling oscillatory relative inflow conditions are used 
to generate the imposed aerodynamic loads on turbosystem 
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4. Transonic aerodynamics 

- For both oscillatory blade motion and oscillatory 
relative inflow conditions, in the absence of 
corresponding two-dimensional transonic aerodynamic 
theories, the transonic blade loads are obtained by 
interpolating between subsonic and supersonic loads. 

For transonic relative Mach numbers at blade tip, the 
subsonic loads are extrapolated. 

5. Causes of aerodynamic excitation 

- These are as described in Section 2.1, Problem 
Description. 

6. Aeroelastic stability at a given operating condition 

- To ascertain that the dynamic response calculations are 
carried out at aeroelastically stable operating 
conditions , the modal flutter analysis capability of the 
present Bladed Disks Computer Program for advanced 
turbopropellers (Refs. 6 and 7) is used. 

2.3 COORDINATE SYSTEMS 


In order to conveniently pose and solve the aerodynamical ly 
forced vibration problem of turbosystems , a number of coordinate 
systems have been defined. Figure 2.3 illustrates these 





coordinate systems for an advanced turbopropeller with its axis 
of rotation mounted at an angle with respect to the tunnel mean 
flow. 

Each of these coordinate systems is described as follows: 

- Tunnel coordinate system 

* This is defined to conveniently specify the velocity 
components of the spatially non-uniform tunnel free 
stream. It can be suitably oriented based on the 
available tunnel data. In the special case of 
aerodynamic excitation in uniform inflow, the tunnel 
coordinate system is oriented such that the X^Z^ plane 
is parallel to the X^Z^ plane of the inertial 
coordinate system as shown in Figure 2.4. The origin of 
the X_^y^Z^ system is arbitrarily located. The 
inclination angle of the turbosystem axis of rotation 
with respect to the tunnel flow also lies in a plane 
parallel to X^Z^ plane . The uniform flow is directed 
along +X^ axis . 

~ XjYjZ^ Inertial coordinate system 

* In the present problem, this coordinate system is used 
to relate the quantities in the tunnel and the basic 
coordinate systems . The orientation of this coordinate 
system is completely arbitrary except for the X^ axis 
to be parallel to , and in the direction of , X^ axis of 
the basic coordinate system described next . The zero 
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reference for time/phase measurements is defined when 
the inertial and the basic coordinate systems are 
parallel, 

All of the following NASTRAN coordinate systems are fixed to 
the rotating turbosystem. 

~ Basic coordinate system 

* This coordinate system has its axis coincident with 
the turbosystem axis of rotation, and directed 
sftward. Location of the origin is arbitrary. The 
XgZg plane contains (approximately) the maximum 
planform of the modelled blade . The definition of this 
coordinate system is consistent with the theoretical 
developments of the 2— d cascade unsteady aerodynamics 
presently incorporated in the Bladed Disks Computer 
Program (Ref. 6) . 

“ (Blade) shank-fixed coordinate system 

* principal advantage of this shank-fixed coordinate 

system is in modelling changes in the blade setting 
angles by a simple 3x3 transformation matrix relating 
to the basic coordinate system. coincides with the 
^isde shank axis . The definition of the coordinate 
system otherwise is arbitrary. 

point location and displacement coordinate 

systems 
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* Any number of such rectangular, cylindrical, or 
spherical coordinate systems can be completely 
arbitrarily defined to locate grid points of the 
NASTRAN model, as well as reguest output at these grid 
points. All of the XgjYgjZg. coordinate systems used for 
output requests collectively form the NASTRAN global 
coordinates system. 

~ Internally generated coordinate system on swept 

chord s 

* This coordinate system is generated within the present 
Bladed Disks Computer Program, and is used to define 
flow and motion properties for the unsteady aerodynamic 
theories on a given swept chord s. It is located at 
the blade leading edge with the x_ directed aftward 
along the chord s. y_ is defined normal to the blade 
local mean surface . 

2.4 DEGREES OF FREEDOM 

The total translational and rotational displacements at any 
fixed point of the rotating and vibrating turobsystem, expressed 
in body-fixed coordinate systems, consist of 

- steady state components due to the steady airloads and 
centrifugal loads , and 

- vibratory components due to the vibratory excitation. 
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superposed on the steady displacements. 


The aerodynamic vibratory response problem of the turbosystem 
is posed herein in terms of the vibratory components of total 
displacements . 


2*5 EQUATIONS OF MOTION AND THEIR SOLUTION 

For an N-bladed tuned turbosystem, with structural coupling 
between blades via a relatively flexible hub, the equations of 


forced motion can be written as (Ref. 4) 

[ ^ } +■ [ [ B" ] + £u"j. 




and 




{p } 




(I) 


I-}., 

Side 


= } 


side 1 


h = 1 , 2 


N 


( 2 .; 


The forcing term on the right hand side of equation (1) is 
entirely due to aerodynamic excitation. Cyclic sector numbers 

and their sides referred to in equation (2) are illustrated in 
Figure 2.5 . 

In seeking solutions for the vibratory displacements u” 
for all n , based on the qualitative and quantitative nature of 
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NOTES 

^ is the angular velocity of the XgY^Zg 

coordinate system w.r.t. the coord . system 

2. Modelled sector is always n=l st. sector . 

3. Sector , and side numbers within a sector , 
increase in the direction of j -o_-t j . 


Figure 2 . 5 Cyclic Sector and Side Numbering 
Convention 



the right hand side forcing functions, the following steps, 
applicable to tuned cyclic structures with rotational cyclic 
symmetry, are considered. 

1. The displacements uP (and the loads ) can be written 

as (Ref. 4) 



This is a standing wave representation wherein the 
coefficients and are, in general, 

functions of time or frequency. For a given circumferential 
harmonic index k, by defining appropriate relations between 
U. and the above equation can be transformed to a 

travelling wave form ( Appendix C ) . The constants kj_ and 
a, are given by. 



QVld 
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2. For a given circumferential harmonic index 
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3. With the use of inter-segment compatibility constraint 
conditions (equation 2) 


{u. } = [ («) 


(7) 


-K 


where u. is an independent displacement vector consisting 
of u. and w. degrees of freedom from the interior and 
side 1 of a cyclic sector. The transformations ^ck. 

functions of the circumferential harmonic k, and 
express the side 2 degrees of freedom in terms of those on 
side 1 . 


4. Equation (5) can then be written as 


n fe 

{^' } = j c6s(«::T k«.) 


} ( 2 ) 


5. The real eigenvalue problem, for a given k, can then be 
stated from equation (1) as 


— 60 


, r —K 

[M j 4- [ K J 




= 0 


(9; 
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6 . 


The eigenvectors ^ J obtained from the solution of 
equation (9 ) , can be used to introduce the modal coordinates 





7. The modal equations of forced motion of the turbosystemr for 
a given circumferential harmonic index k , can then be written 
from equation (1) 

[M ]{t ) *■ [ ® ^ 

where 

3 (rs) cov^ta. 


=r K K 
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, K T Y. 
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C^cfeT CQ [G^jl Cq^DL* 


^ Oj) 


j 

—— 

8. For a given circumferential harmonic index k, P and 
— fe.5 

P are the circumferential harmonic components of the 
total external excitation. Such excitation due to 
aerodynamic sources is discussed in Section 3. 


9. The generalized oscillatory aerodynamic reaction matrix 
Q can be written as 

p 


where 




(\S') coy^-td. 
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Os) 

are the 'cosine' and 'sine' component mode shapes of 

- K 

eigenvectors 0 Equation (14) is rewritten as 


For turbosystem structures with flexible hub/disk/ 


a) Q and Q , . exist when the circumferential harmonic 

index k 0 and qfc N/2 when N, the total number of cyclic 
segments in the structure, is even, and 

c 

b) only O.. exists when k = 0 or N/2 , N even . 

For turbosystem structures with rigid hub/disk , each cyclic 

segment of the structure behaves structurally independent of 

its adjacent segments . Degrees of freedom at segment 

boundaries are completely constrained to zero . The only 

possible structural modes are those akin to k = 0 modes with 

c. 

fixed inter-segment boundaries . Only Q . . exists . 


Derivation of < 2 . . is discussed in detail in Appendix A. 
Supersonic chords are considered as examples . The derivation 


of Q. . for subsonic chords alone was discussed in Ref . 

C* S 

Q . . and Q, . . are generated , in turn , by considering 
— ks 

and ^ , in turn, to represent the structural modes in 


6 . 

— kc 
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equations (A29) and (A31) of Appendix A. 


10. Equation (12) can now be solved for \ . Substitution in 

M Is. 

equation (11) / and equation (8) yields u ' . 

11. Repeating steps 2 through 10 for all applicable 
circumferential harmonic indices, and substitution in 

■A 

equation (3) , result in for all n . 

12. Other dynamic responses such as stresses, etc., can be 
obtained for all sectors of the turbosystem by current 
NASTRAN procedures. 

The procedure described above for the solution of 
circumferential harmonic components of dynamic response, for a 
given circumferential harmonic index , has been implemented in the 
April 1984 release of NASTRAN on the CRAY 1-S computer system at 
NASA LeRC (Ref. 9) . 

An overall flowchart of the solution procedure is shown in 
Figure 2.6 . 
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3. APPLIED OSCILLATORY AIRLOADS 


3.1 GENERAL 

Generation of applied oscillatory aerodynamic loads on the 
blades of the turbosystems due to oscillatory relative inflow is 
briefly discussed in this section. Complete details of the 
airloads generation process are reported in Ref. 8. 

A rotating turbosystem placed in a spatially non-uniform 
inflow, or in a uniform inflow with its axis of rotation 
misaligned with inflow direction , experiences oscillatory 
relative inflow. 

The resultant oscillatory airloads are computed by the use of 
subsonic and supersonic 2-d cascade unsteady aerodynamic theories 
(Refs . 10 and 11) . Accordingly, the aerodynamic model of the 
blade comprises non-intersecting strips as shown in Figure 2.2. 
The swept blade of an advanced turboprop is shown as an example . 
The analytical development is equally applicable to unswept 
blades . 

The computation of the oscillatory airloads is carried out in 
the following four steps: 

1 . For a given chord , sinusoidal gust amplitudes , frequencies , 
and other aerodynamic excitation parameters ( Appendix C ) , 
are determined based on the relative inflow variations as the 
blades go through one revolution . 
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2. Subsonic or supersonic theory with appropriate gust inputs 
from step 1 is used to determine oscillatory pressure 
distributions on the chord. 

3. The unsteady pressure distribution on the blade strip 
associated with the chord is transformed to loads at the 
finite element structural grid on the chord. 

4. Steps 1 through 3 are repeated for all chords spanning the 
blade . 

After establishing the relative inflow velocity to be either 
subsonic, transonic, or supersonic from step 1, step 2 is merely 
an application of the appropriate theory. In case of transonic 
relative inflow, the airloads for that chord are interpolated 
from adjacent non-transonic chords . Steps 1 and 3 are discussed 
further in the following sections . 

3.2 GUST AMPLITUDES AND FREQUENCIES OF OSCILLATORY 
RELATIVE INFLOW 


The problem of determining the amplitude and frequency 
contents of oscillatory relative inflow gusts can be accomplished 
in four steps : 

1 . At the leading edge point A of any chord 5 (Figure 3.1), 
the relative inflow velocity is determined as a function of 
time t (at discrete time instances) during one revolution of 
the turboprop rotating at a constant angular velocity si.. 
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2. At each time step, the relative inflow velocity can be 

resolved into ^ 5- r and Zy components of the local 

^ <s 

chord coordinate system (Figure 3.1) . 

- Component along is the cascade relative inflow 
velocity. Component along 3^ is the gust velocity. 
Component along "Zy is the radial velocity, which is 
ignored in 2-d cascade aerodynamics. 

3. A Fourier series decomposition in -a-t of 

a) the velocity component along x„ yields the cascade 
mean relative inflow velocity as the constant term of 
the Fourier series (along AB, Figure 3.1). This is used 
to determine the relative Mach number and select the 
appropriate aerodynamic theory. 

b) the velocity component along yields the necessary 
sinusoidal gust amplitudes and frequencies which are 
subsequently input to appropriate oscillatory inflow 
aerodynamic theories . 

4. Steps 1 through 3 are repeated for all chords spanning the 
blade to determine the oscillatory pressure distribution over 
the entire blade surface . 

Referring to Figure 3.1, the relative inflow velocity at the 
leading edge point A on the chord s , expressed in the chord 
coordinate system, is 
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equation (1) reduces to 
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For clarity of presentation/ further development is separated 
into that for uniform and non-uniform inflows. 


3.2.1 Uniform Inflow 

In the uniform inflow case, the turbosystem axis of rotation 
is inclined at an angle with the uniform inflow velocity as 

shown in Figure 2.4. The transformation from the tunnel to the 
inertial coordinate system (equation (4) ) becomes 


r ^ "1 


L 


cosY 

o 

si^Y 


-Si^Y 

o 

coiY 


(9) 


with the uniform inflow velocity given by 


To»^v\e \ 

{wa} = ^ 


(wa). 


o 

D 


(io) 
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Substitution of equations (10), (9), (5) and (3) in equation 
(7) yields the constant (independent of -a-t ) part of (VA) — , 
and the oscillating (function of ) part of (VA) - as 


co«(va)_ = p 
X. ' 11 




os£ = (u/a)^ si‘«y sm -n.t +• p cex-f^'t'J 

The equivalent complex representation of equation (12) 
CVA)_ = (WA)^^ e" 

Witt, toA s: - 


O’-; 


is 


0?; 

(lA) 


Equation (11) defines the cascade relative inflow velocity for 
the chord s , whereas equation (13) describes the incoming gust 
velocity at the frequency _n. , 


3«2.2 Non-Uniform Inflow 

Equation (7) , in its most general form, yields the expressions 
for the cascade relative inflow velocity and the gust velocity. 
The cascade inflow velocity is 


(^A)_ = (Wa)^^ 4. 9,3 (Wa)^_^ 

^ A (O A) n 

s ' ^ > 


(15^^ 


where g^j 's are the elements of the transformation G (equation 
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Both (VA)_ (equation (15) ) and (VA) „ (equation (19) ) are 
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periodic functions of the azimuthal angle -il-fc , with a maximum 
period of lir as the blade completes one revolution in spatially 
non-uniform absolute inflow field. Each of these velocities can 
be expanded as finite Fourier series 

P 

{ C^t) = 2 

The coefficients f^^ , determined by knowing 

f at M equally spaced intervals between 0 and itt . The upper 
limit of harmonics , P , is given by 



P 


( M - O / 2 H o4^ 

/ '' 

H /2- ^ M evCA , 




In the non-uniform inflow case, the f ^ component of (VA)^ 
is taken as the constant part of the cascade relative inflow 
velocity. For a selected harmonic p , the oscillating part of 
(VA) y is written as 


osc 





with 





/ % 


i2e) 


The excitation frequency correspondingly is hJX. 
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3.3 GRID POINT LOADS FROM OSCILLATORY PRESSURE 

distribution along chord 

The subsonic and supersonic aerodynamic routines compute the 
oscillatory pressure distribution at a number of points between 
the leading and trailing edges of a given chord. These points 
are generally distinct from the structural grid points at which 
the applied airloads are desired. 

To obtain the loads at structural grid points, the chord is 
divided into a number of segments as shown in Figure 3.2. For 
each of the segments, the loads at its ends, directed along 
+ ^ , are calculated. As an example, for the G^ segment. 



with the given chord . 


Furthermore, the load at grid point G^, for instance, 
comprises the algebraic sum of those contributed by segments 
G ^ G 2 _ and G * 

The applied oscillatory airloads generation procedure has been 
implemented in a stand-alone (independent of NASTRAN system) 
computer program, AIRLOADS, which is described in detail in Ref. 

8 . 
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4. APPLICATION EXAMPLES 


4.1 GENERAL 

The theoretical development of the preceding sections, as 
implemented in the Bladed Disks Computer Program in the April 
1984 release of NASTRAN on the CRAY 1-S computer system at NASA 
LeRC, was used to conduct a total of eight examples on the SR-3 

and SR-5 advanced turboprops. These examples were selected by 
NASA. 

4.2 EXAMPLES 

The SR-3 and SR— 5 examples are listed by operating conditions 
in Table 4.1. These conditions were selected by NASA to cover a 
range of wind tunnel test conditions and provide a degree of 
commonality between examples to cut down on the computational 
steps. 

At each of these operating conditions, it was desired to 

1. calculate the natural frequencies and mode shapes of the 
turboprops, 

2. examine the aeroelastic stability, and 

3. determine the one-per-rev stress response at the strain gage 
locations shown in Figures 4.1 and 4.2 for the SR-3 and SR-5 
turboprops, respectively. 
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Table 4 . 1 Operating Conditions for SR- 3 and SR-5 Examples 


Turboprop 

Blade 

NASA 

Test 

Rdg. 

No. 

Freestream 

Blade 

Setting 

Angle 

P 3/4 (deg) 

Prop, axis 
tilt 
(deg) 

Prop. 

RPM 

Mach 

No. 

Velocity 
(ft/ sec) 

Density 

X 10 , 

(Ibf-sec / 
ft^) 

SR- 3 

190 

.353 

398 


48.9 

8 

8000 

273 

.798 

873 

mm 

60.8 

2 

8000 

277 

.795 

871 

1.9053 

60.8 

4 

7000 


.795 

871 

1.9060 



60.8 

4 

8000 

SR-5 

8508 

.801 

885 

n 

70.8 

2 

5500 

8511 

.797 

881 

1.8313 

70.8 

5 

5500 

8607 

.360 

411 

2.3057 

60.8 

1 

12 

6000 

8610 

.361 

411 

2.3062 

60.8 

3 

6000 






















































4.3 STRUCTURAL AND AERODYNAMIC MODELS 


NASA provided the NASTRAN finite element structural and 
aerodynamic models of the SR-3 and SR-5 turboprops shown in 
Figures 4.3 and 4.4/ respectively. The aerodynamic models were 
defined as subsets of the structural grid. 

Both models, with homogeneous material properties of Young's 
modulus = 1.6 X 10 psi , Poisson ' s ratio = .35 , and mass 
density = .16 Ibm/ in , were built using the TRIA2 triangular 
plate element with bending and membrane capabilities. 


4.4 BOUNDARY CONDITIONS 

For the purposes of the present analyses , the turboprop hub 
was considered to be rigid as compared to the blade flexibility 
for both SR-3 and SR-5. This mechanical decoupling between blades 
was achieved by completely constraining the degrees of freedom at 
the bottom of the blade shank, and by not modelling the hub . 


4.5 APPLIED OSCILLATORY AIRLOADS 

The applied airloads were generated by using the AIRLOADS 
program discussed in Section 3 and Ref. 8. This program is 
capable of reading the NASTRAN structural and aerodynamic grid 
data, and the operating conditions to generate the oscillatory 
airloads . Upon user 's request , the program also produces two 
output files for the oscillatory applied airloads for convenient 
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Figure 4.3 NASTRAN Structural and 

Aerodynamic Models of SR- 


S-HS JO 
pup 


spppow OTuiPuApouav 
XPjnaonu:^S MVHXSVM 


aunbx^ 





inclusion into NASTRAN Case Control and Bulk Data Decks. The 
streamline geometry and flow data generated by the AIRLOADS 
program was also useful in conducting aeroelastic stability 
analysis. 

4.6 ANALYSIS STEPS FOR EACH OPERATING CONDITION 

A series of procedural steps was setup to systematically 
conduct all of the analysis required for the turboprops at each 
of the operating conditions listed in Table 4.1. An overall 
flowchart of these steps is shown in Figure 4.5. 

The sequence of four runs shown is arbitrary to an extent , 
and could be rearranged as per user 's needs . 

Run 1 . Differential/St iff ness analysis was conducted at the 
i^pm with the blade set at the selected blade setting 
angle. Displacement Approach RF 4 was used , and the total 
{elastic plus differential) stiffness was saved as K TOTAL . 

Run 2/3. Aeroelastic stability of the turboprop at the 
selected operating point was determined using the Aero Approach 
RF 9. V-g , V-f plots were used to determine stability. 

If desired / the AERO and STREAMLi bulk data cards required for 
this run can be obtained from the AIRLOADS program. Run 3/2. 

Natural frequencies and mode shape plots can be optionally 
obtained in this run. 

Run 3/2 . At the given operating condition with the prop axis 
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Figure 4.5 Overall Flowchart of Analyses Steps- to Conduct 
Each of the SR— 3 and SR— 5 Examples 










set at an angle with the uniform inflow, the applied oscillatory 
airloads distribution on the modelled turboprop blade was 
generated using the newly developed AIRLOADS program. At user's 
option, individual output files for direct inclusion in the Case 
Control and Bulk Data Decks of the subsequent response Run 4 were 
also created. 

Run 4 . Finally, modal forced vibration analysis was conducted 
using the newly developed Displacement Approach RF 8 alter 
package, MFVAAET (Modal Forced Vibration Analysis of 
Aerodynamically Excited Turbosystems) . 

Element stresses in element coordinate systems for elements 
neighboring a strain gage location were transformed along strain 
gage orientation as follows; 


r S — S 

= 1 + \ — 

S — S’ 

C — — 1 si/v 18 — S cos 1^ 

sKea-1* 2. / 


where 9 is the angle between the strain gage axis and the 
x~axis of the element coordinate system. 
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5. RESULTS AND DISCUSSION 


5.1 GENERAL 

Results of the SR-3 and SR-5 examples described in the 
previous section are presented in Table 5.1 (SR-3 and SR-5), 
Figures 5.1 through 5.15 (SR-3), and Figures 5.16 through 5.32 
(SR-5) . 

For all of the contour plots shown in these Figures, the 
contours for the translational displacements along the Basic 
coordinate system Y axis are plotted within the UNDEFORMED shape 
outline. The outline for the DEFORMED shape accounting for the 
total translational displacements is also shown. For all mode 
shape plots, the contours represent zero displacement along the 
Basic Y axis . 

Following discussions are grouped by analyses . 


5.2 DIFFERENTIAL STIFFNESS 

Results of differential stiffness analysis of SR-3 test 
reading no. 190, and SR-5 test reading no. 8508 are shown in 
Figures 5.1 - 5.2, and 5.16 - 5.17 respectively. 

Convergence of the non-linear part of the analysis was 
achieved after two iterations for all operating conditions . 

In the present context , the differential stiffness analysis 
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Table 5 . 1 Natural Frequencies of SR- 3 and SR-5 Turboprops 


Turboprop 

Blade 

NASA 

Test 

Rdg. 

No. 


Prop . 
RPM 

Natural Frequencies, Hz. 


Mode 2 

Mode 3 

Mode 4 


SR-3 

190 

48.9 

8000 

227 

431 

673 

759 

915 

273 

60.8 

8000 

220 

423 

666 

757 

911 

277 

60.8 

7000 

210 

417 

652 

743 

905 

278 

60.8 

8000 

220 

423 

666 

757 

911 

SR-5 

8508 

70.8 

5500 

140 

275 

558 

650 

718 

8511 

70.8 

5500 

140 

(166) * 

275 

(287) 

558 

650 

718 

8607 

60.8 

6000 

150 

281 

564 

658 

723 

8610 

60.8 

6000 

150 

(173) 

281 

(288) 

564 

658 

723 


* 


Experimental Values 













MAX-CEF. = 0. 


36: 



SR3 RESPONSE TO 1 PER REV OSC. AIR 
NASA TEST READING NO. 190 
DIFF. STIFF. ANAL.— LINEAR SOLUTIOf 
STATIC DEFOR. SUBCASE 1 LOAD 1 


Figure 5 . 1 



SR3 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NO. 190 
OIFF. STIFF. ANAL.— NONLINEAR SOLUTION 
STATIC DEFOR. SUBCASE 2 LOAD 1 

Figure 5 . 2 SR- 3 Dif 

Analysis 
Solution 
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MAX-OEF. = 0.55044185 



SR3 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NO. 190 
K - 0 MOOES- OSCILLATORY AIRLOADS PRESENT 
MODAL DEFOR. SUBCASE 1 MODE 1 FREQ. 227.1975 

Figure 5.3 SR-3 Mode 1 
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SR3 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NO. 190 
K » 0 MODES- OSCILLATORY AIRLOADS PRESENT 
MODAL OEFOR. SLBCASE 1 MODE 2 FREQ. 431.7661 

Figure 5.4 SR-3 Mode 2 
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MAX-DEF. = 0.27590971 



SR3 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NQ. 190 
K - 0 MOOES- OSCILLATORY AIRLOADS PRESENT 
MODAL DEFOR. SLBCASE 1 MODE 3 FREQ. 673.0511 


Figure 5 . 5 SR- 3 Mode 3 
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SR3 RESPONSE TO 1 PER REV OSC. Aira_OADS 
NASA TEST READING NO. 190 
K = 0 MOOES- OSCILLATORY AIRLOADS PRESENT 
MODAL DEFOR. SLBCASE 1 MODE 5 FREQ. 915.2528 


Figure 5 . 7 
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Figure 5 . 8 


SR-3 V-g, V-f Curves 
— Test Reading No. 190 
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Figure 5.9 SR- 3 V-g, V-f Curves 

— Test Reading No. 273 
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Figure 5 . 10 SR- 3 V-g, V-f Curves 

--Test Reading No . 277 
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Figure 5-11 


SR-3 V-g, V-f Curves 
— Test Reading No . 278 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 


SR3: NASA TEST READING NO. 190 
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Figure 5.12 SR- 3 One-Per-Rev Stress 

Comparison — Test Reading 
No. 190 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 
SR3: NASA TEST READING NO. 273 


X TEST VALUE 


RANGE COVERED BY FINITE 
ELEMENTS NEIGHBORING THE 
STRAIN GAGE LOCATION 
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5.13 SR— 3 On0— Per— Rev Stress 
Comparison — Test Reading 
No. 273 
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MAGNITUDE OF OSCILUTORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 
SR3: NASA TEST READING NO. 277 
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Figure 5.14 SR- 3 One-Per-Rev Stress 
Comparison — Test Reading 
No. 277 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 
SR3: NASA TEST READING NO. 278 
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Figure 5.15 SR- 3 One-Per-Rev Stress 

Comparison — Test Reading 
No. 278 
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SR5 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NO. 8508 
K ~ 0 MOOES- OSCILLATORY AIRLOADS PRESENT 
MODAL DEFOR. SUBCASE 1 MODE 1 FREQ. 140.09^18 


Figure 5.18 SR-5 Mode 1 
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MAX-OEF. = 1.00000000 



SR5 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NO. 8508 
K - 0 NODES- OSCILLATORY AIRLOADS PRESENT 
NODAL DEFOR. SUBCASE 1 NODE 2 FREQ. 275.8804 


Figure 5.19 SR-5 Mode 2 
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SR5 RESPONSE TO 1 PER REV OSC. Aira_OAOS 
NASA TEST READING NO. 0508 
K - 0 MOOES- OSCILLATORY AIRLOADS PRESENT 
MODAL DEFOR. SUBCASE 1 MODE 3 FREQ. 558.1463 


Figure 5.20 SR- 5 Mode 3 
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SR5 RESPONSE TO 1 PER REV OSC. AIRLOADS 
NASA TEST READING NO. 8508 
K - 0 MOOES- OSCILLATORY AIRLOADS PRESENT 
MODAL DEFOR. SUBCASE 1 MODE 4 FREQ. 660.2884 


Figure 5.21 SR- 5 Mode 4 
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SR5 RESPONSE TO 1 PER REV OSC. AIRLOADS 

NASA TEST READING NO. 0500 
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MODAL DEFOR. SLBCASE 1 MODE 6 FREQ. 718.7901 
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Figure 5.23 SR- 5 V-g, v-f Curves 

— Test Reading No. 8508 


70 







O^IGIJVAL 

OF POOR 


13 

quality 


e.B l.a 1.6 2.8 2.B 3.8 3.6 4.8 4.6 E4 



Figure 5.24 SR-5 V-g, V-f Curves 

— Test Reading No. 8511 
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Figure 5.25 SR-5 V-g, V-f Curves 

— Test Reading No. 8607 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 
SR5: NASA TEST READING NO. 8508 
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Figure 5 . 29 SR-5 One-Per-Rev Stress 
Comparison — Test Reading 
No. 8508 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 
SR5; NASA TEST READING NO. 8511 
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Figure 5.30 SR- 5 One-Per-Rev Stress 
Comparison — Test Reading 
No . 8511 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 


SR5: NASA TEST READING NO. 8607 
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Figure 5.31 SR- 5 One-Per-Rev Stress 
Comparison- -Test Reading 
No. 8607 
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MAGNITUDE OF OSCILLATORY STRESS (PSI) 


COMPARISON OF ANALYTICAL AND TEST STRESSES 


SR5: NASA TEST READING NO. 8610 
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Figure 5.32 SR- 5 One-Per-Rev Stress 

Comparison- -Test Reading 
No. 8610 
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was conducted principally to save the total (elastic plus 
differential) stiffness matrix for use in subsequent aeroelastic 
stability and dynamic response analyses. 

5.3 NATURAL FREQUENCIES AND MODE SHAPES 

The first five natural frequencies for all SR-3 and SR-5 
operating conditions are listed in Table 5.1. These frequencies 
were obtained by including the effects of differential stiffness 
and centripetal softening terms. 

The only available experimental frequencies are also shown in 
the Table . 

The first five mode shapes for SR-3 test reading no. 190, and 
SR-5 test reading no. 8508 are shown in Figures 5.3 - 5.7, and 
5.18 - 5.22 , respectively. These mode shapes were determined at 
the operating speeds of the turboprops. Inserts of holograms 
from the bench modal tests (zero rpm) are also included in these 
Figures . Although the nodal contour patterns can not be directly 
(and exactly) compared, due to the rpm differences and the Basic 
coordinate system Y axis not being normal to the blade surface at 
all points , a broad-based correlation between the calculated and 
observed mode shapes is indicated in these Figures . 


5.4 AEROELASTIC STABILITY 

Velocity-damping (V-g) and Velocity-frequency (V-f) curves 
summarizing the complex eigenvalues of aeroelastic stability 
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analysis for each of the four operating conditions of the SR-3 
turboprop are shown in Figures 5.8 through 5.11. 

The first five natural inodes were included in the stability 
analysis conducted by the KE method in NASTRAN. The reduced 
frequencies were based on the natural frequencies of these modes, 
and the semichord and cascade relative inflow velocity of the 
reference chord (Figure 4.3). No structural damping was assumed. 
The cascade relative inflow velocity on the reference chord 
defining the operating point is shown by a vertical line in each 
of the V-g, V-f plots. The aeroelastic stability of the 
operating point is established by a negative damping value for 
each of the structural modes at the operating point velocity 
(vertical line) . 

This is observed to be true in all of the SR-3 operating 
conditions (Figures 5.8 through 5.11) . 

The V-g, V-f curves for the four SR-5 operating conditions are 
presented in Figures 5.23 through 5.26. 

In these cases too, the first five structural modes were 
considered in the aeroelastic stability calculations . No 
structural damping was assumed . 

While the operating points 8607 and 8610 (Figures 5.25 and 
5.26) are seen to be stable , aeroelastic instability is indicated 
at the 8508 and 8511 test conditions (Figures 5.23 and 5.24) by 
the positive damping value in structural mode 1. To highlight the 
damping distribution in this mode near the operating velocity. 
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these two cases were further analyzed by selecting a finer 
reduced frequency range (.1 to .3) about the first structural 
mode frequency. The resulting V-g, V-f plots are shown in 
Figures 5.27 (8508) and 5.28 (8511) . A positive aerodynamic 
damping value of about 2% is indicated in mode 1 for both of 
these curves. With the assumption of 2-3% structural damping 
present in the turboprop blades, the 8508 and 8511 test 
conditions can be treated as marginally subcritical operating 
points. A uniform structural damping value of 3% was used in the 
response calculations at these two operating conditions. 

5.5 VIBRATORY RESPONSE TO AERODYNAMIC EXCITATION 

The aerodynamic excitation in all of the eight examples is due 
to the turboprop axis of rotation being misaligned with the 
uniform inflow. Hence, the only excitation frequency is the 
one-per-rev value for each of the operating conditions . 

Magnitudes of the one-per-rev oscillatory stress response of the 
SR- 3 and SR- 5 turboprops at the strain gage locations with 
available wind tunnel data are shown in Figures 5 . 12 - 5.15, and 
5.29 - 5.32 respectively. (Strain gages are identified in 
Figures 4.1 (SR-3) and 4.2 (SR-5) ) . 

A range of stress values covered by the finite elements 
neighboring a strain gage location is plotted for direct 
comparison with the measured value . 

The first five structural modes in each case were included in 
the response calculations . 
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In regards to the comparison of calculated and measured 

stresses, the following points are noted; 

1. In modal forced response analysis, the calculated stress 
response is essentially determined by the following factors; 

a. modal frequencies, 

b. modal stresses, 

c. modal damping, 

d. excitation frequencies, and 

e. excitation levels. 

2. If the excitation frequencies are well separated from the 
modal frequencies , the modal damping .exerts insignificant 
influence on the stress response . 

3 . If only one excitation frequency exists, and is placed well 
below the first modal frequency, the principal , and 
significant , contribution to the response stresses comes from 
the modal stresses of the first mode . 

4 . The contributing modal stresses are 'scaled ' by the 
excitation levels to the calculated response stresses . 

5. Keeping in view the above 4 points , the following 
observations are made in the eight SR- 3 and SR- 5 examples; 

a. The only one-per-rev excitation frequency in each 

case is well below the first modal frequency as 
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summarized in Table 5.2 . 


b. This implies the insignificance of modal damping, 
and the significance of modal stresses of the first 
structural mode toward computing the response stresses . 

Figure 5.33 illustrates the striking similarities 
between the calculated stress responses of all four SR-3 
cases . Gage-to-gage stress distribution for each case 
(Figures 5.12 - 5.15) is normalized with the mean stress 
of gage 1 in that case. These stress distributions are 
also noted to closely resemble those exhibited by the 
modal stresses in the first structural mode . 

Figure 5.34 shows a similar relation between the 
calculated stress responses of all four of the SR- 5 
cases . 

c. The one-per-rev applied oscillatory airloads 
distribution effectively scales the modal stresses in 
the first mode to the calculated response stresses . 

While an across the board inspection of comparisons of 
calculated and measured stresses in all eight cases reveals a 
generally good correlation, some specific deviations are 
noted as follows : 

a. In all SR— 3 cases , gage 3 measurements are higher 
than predictions . This may be due to discrepancies in 
strain measurements or strain-gage orientation . The 
finite element model (Figure 4.3) does not reflect 



Table 5 . 2 Comparison of Excitation and Modal Frequencies 


Turboprop 

Blade 

NASA 
Test 
Rdg . 
No. 

Excitation 

Frequency^ 

Hz, 


Mode 1 

Mode 2 

SR- 3 

190 

133 

227 

431 

273 

133 



277 

117 

210 

417 

278 

133 

220 

423 

SR-5 

8508 

92 


275 

8511 

92 

140 

275 


100 

150 

281 


EM 

100 

150 

281 
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Figure 5.33 SR-3 One-Per-Rev Calculated 

Stress Response Comparison Between 
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SR-5: ALL CASES 


1-00 T 


X 4- - o 


w 

D 


> 


.804 


SR”5 Test Reading 


X 

8508 


8511 


8607 

o 

8610 


w 

o 

o 

>8 

Q 

W 


£S3 

H 




o 


.60 4 


.40 1 


.201 


0 


^ ^ ^ o 


GAGE 1 GAGE 3 GAGE 5 GAGE 6 


STRAIN GAGES 


Figure 5.34 SR- 5 One-Per-Rev Calculated 

Stress Response Comparison Between 
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abnormal element connections or boundaries in this 
region . 

b. Similar observations can be made for gage 6 in SR-5 
cases . 

c. Test cases with a combination of low (.35 - .36) 

Mach numbers (SR-3 190, SR-5 8607 and 8610), high tilt 
angles (SR-3 190 and SR-5 8607) , and very highly swept 
blades (SR-5 8607) may be more sensitive to the 
aerodynamic modelling of the tip regions of the blade. 
This, in turn, affects the applied oscillatory airloads, 
ultimately altering the vibratory stress response . 
Spacing, number, and orientation of streamlines in such 
regions of the blade are some of the variables of the 
aerodynamic model which will provide a good basis for 
any sensitivity investigations . 
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6. CONCLUSIONS 


1. Theoretical developments of a new capability in NASTRAN for 
modal forced vibration analysis of aerodynamically excited 
turbosystems have been described. 

2. The capability has been successfully applied to eight NASA 
wind tunnel test conditions of the SR-3 and SR-5 advanced 
turboprops for predicting the vibratory response of the 
turboprop blades. 

3. NASTRAN pre-processing capability of the stand-alone AIRLOADS 
program has been successfully utilized in generating the 
applied oscillatory airloads on the swept blades of the 
advanced turboprops. 
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7 . RECOMMENDATIONS 


The application examples conducted in the present work have 
been for single-rotation advanced turboprops at a number of 
operating conditions. Two or three of these eight cases have 
indicated a need for analyzing refined aerodynamic models in 
improving predictions. 

It would prove useful tOr 

1. study the aerodynamic model variations as they affect the 
applied airloads and response predictions, and 

2. apply the capability to analyze forced vibrations of other 
turbosystems such as counter-rotating turboprops and 
axial-flow turbomachines. 
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APPENDIX A 


SWEEP EFFECTS IN SUPERSONIC AIRLOADS DUE TO 
OSCILLATORY BLADE MOTIONS 


A.l GENERAL 

The generalized modal aerodynamic force matrix 0.^1 due to 
oscillatory motions of the swept blades of advanced 
turbopropellers is derived. The two-dimensional supersonic 
cascade unsteady aerodynamics program of Ref. 11 has been 
modified to include the effects of blade sweep. The 
modifications are similar to those carried out for subsonic flow 
{Ref. 6) . The blade is spanned by a number of non-intersecting 
chords selected normal to any spanwise reference curve such as 
the blade leading edge (Figure A.l) . The modified subsonic or 
supersonic two-dimensional cascade theory is applied on each of 
these chords to determine the generalized aerodynamic forces 
acting on the associated strip. The strip results are added to 
obtain the blade aerodynamic matrix . This Appendix discusses the 
sweep changes in supersonic aerodynamics . Those for subsonic 
flow are contained in Ref. 6 . 

In order to conform to the existing computational scheme in 
the Bladed Disks Computer Program, and thus reduce extensive 
recoding , the chordwise generalized aerodynamic matrices are 
first compu ted for predefined , simple , chordwise aerodynamic 
modes (Figure A. 2) . The chordwise structural modes are then 
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introduced via Fourier expansions to finally obtain the blade 
generalized airforce matrix. 

A. 2 DEFINITIONS 

Figure A.l illustrates some of the definitions pertinent to 
incorporating sweep effects in the 2-d cascade programs. A_B_, 
AB and A^B^ represent three successive chords with points A's on 
the leading edge. For the chord AB, at any operating condition 
WA represents the absolute inflow velocity while AO (= ^ x RA) 

is the blade (tangential) velocity. 

In the plane WAU, VA = WA - AD represents the relative inflow 
velocity. 


AI is parallel to the axis of rotation, and does not , in 
general , lie in the WAD plane . 


AD is the projection of AC (BA extended arbitrarily to C) in 
plane lAD. 


The sweep angle A. , and the stagger angle ^ are 
positive as shown . 


A local coordinate system xyz is defined at the leading edge 
point A of the chord AB such that x is directed along AB. y is 
defined normal to the 'mean ' surface containing the points A_ , A, 
A^, B, and B_. The unit vector along y , for the sense of ^ 

shown in Figure A.l is given by 


j 



(a.b^)x (ab) 


(AS) X 

(^) X (^_)j 


(A1) 
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Modal translations along y and rotations about “x are used in 
deriving the generalized airforce matrix (Section A. 6) . For the 
opposite sense of rotation, xyz is defined to be left handed with 
y reversing direction. 

The shaded area about the chord AB represents the strip of 
integration associated with AB. The length of the strip is given 
by AB, the widths at the leading and trailing edges are 
respectively given by 

V = i 

= i t 

A 

where k is the unit vector along z". 

A. 3 BLADE DEGREES OF FREEDOM 

For aerodynamic computations , consider the blade to be spanned 
by S chords with G structural grid points on each chord . The 
blade degrees of freedom £u j- can then be partitioned as 

^1 

s 

where ^ 
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M = 


(A + ) 


u represents the 3 translational and 3 rotational degrees 
of freedom at grid point g on chord s. 


A. 4 NORMAL DISPLACEMENT, DOWNWASH AND PRESSURE DISTRIBUTION IN TERMS 
OF CHORDWISE AERODYNAMIC MODES 


The normal displacement (along + y) at any point x on the 
chord s can be expressed in the local (chord) coordinate system 
as a linear combination of the chordwise aerodynamic modes : 


* .-k } = A 


u ) [ (as) 


where is the semi-chord , 

s 


L$'J = L' £ 


/n-l. .TT. X.\ . / 9-Z .IT. £ \ (6i-l f . 

(A6) 


represents the chordwise aerodynamic mode shapes , and 




{ki) 


is the aerodynamic modal participation vector . 
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The downwash at x can be written as 




where V - is the chordwise cascade relative inflow velocity 
directed along CA. 

Substituting for U from eguation (A5 ) , the downwash can be 
written as 




sv\ 




t- tav\A L ^ J 




L^t 


(a9; 


where 


= 


r 






si 


(au) 


and the reduced frequency 


k_ = 60 1 - / V 


(A11) 


Corresponding to the downwash boundary conditions associated 
with each of the elements of the modal participation vector 

T ^ . . ?2- 

i its derivative ^ the modified 2^-d supersonic 

cascade unsteady aerodynamics program computes the pressures on 
and surfaces at P discrete points distributed along 

the chord s. (The contribution due to the variation of the 
perturbation velocity potential in the z direction is 
neglected.) The distributions of these P points on the y_ and 
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surfaces of the chord are distinct in general. Hence, the 
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where 




whereby the generalized chordwise aerodynamic force matrices are 


With I r;,j = I x;p ^ 

ii;^3 = 


^ 


I (azi ) 


> (Ai^) 


J 


Elements of matrices and have been evaluated as 

chordwise integrals, and are presented in Appendix B. 


A. 6 TRANSFORMATION BETWEEN AERODYNAMIC AND STRUCTURAL MODAL 
COORDINATES 


Equation (A5) can be used to express the normal displacement 
at the G structural grid points on chord s as 

{ "} = 

In terms of structural modal coordinates , w. can also be 
written as 



9 ^ 
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Comparison of equation (A23) with equation (A24) yields 


= ir 


(a 2.6^ 


Differentiation with respect to 2 results in 

,sl 


r9'^ 


T } ‘ 5 • i- H rap • K; j] ty . 


where 


r n r 3- ^ "T* r -N ^ w Uf^J 

= tv tis V ^ 


(Azg^ 


!o£^I 

The modal displacement matrix ^ and its derivative 

^ i !oC€k j 

^ p are derived from the global modal matrix p as 

35 

follows : 


SH JorgL 


^ global 

t 1 - ] 
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(AZ^) 


where each of the i modal columns of jzJ consists of 
the three modal translations at each of the G computing station 
grid points on the chord s . 


sni 

The transformation G is given by 
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0 1 o 

0 1 0 


^ 1 
[r-^l 


• 

L 1 i « - 

^ SX J 




-7 


0 1 O 




where and represent the coordinate transformations from 
the NASTRAN basic system to the local (chord) xyz system and the 
global (displacement) system, respectively. 


Similarly, 


Ilil. 3 = Ev]!^V 


(A3l) 


glcVtixi . 

p consists of 


where each of the i modal columns of 
the three modal rotations at each of the G computing station grid 
points on the chord s. 


The transformation G 


JV5 2. 


is given by 
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GENERALIZED AIRFORCE MATRICES FOR CHORDWISE STRUCTURAL MODES 


Introducing relations (A26) and (A27) in equation (A20 ) , the 
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-T 


[Cl = i- [C] ^rvijh-) 


= i />y,^ I - ro.';i 


si z2 

Combining Q and Q , the generalized airforce matrix for 
chord s is defined as 




(azg) 


A. 8 BLADE GENERALIZED AIRFORCE MATRIX 


The virtual work done by the aerodynamic forces on the blade 
is the sum total of the virtual work done by all the chordwise 
aerodynamic forces on their respective strips . Thus 


.hloide 


= xw 


- ^ C,i= y <SL?- U\-} 

z ‘ i ^ ref ^ L J ^ ^ 


whereby the generalized airforce matrix for the swept blade of 


the advanced turbopropeller can be written as 

= I ^ V-"' = i/-’" , [5 1 

L J r S vef s^ref / J x t .ref L Xl-J 


(AZ8) 
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APPENDIX B 


CHORDWISE WEIGHTING MATRICES (SUPERSONIC) 


The chordwise weighting matrices and appearing in 

equations (A22) of Appendix A are evaluated. 
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The six typical elements of I„t, are as follows, with 
y= 3c/£j : 
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s Oiy\J (f + s) £\/&v\ 
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The six typical elements of are as follows, with y = x/jB.^: 

I. = f = ' 1 . 


b) as given by typical element (a) of I^^j, 
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(r-7.) 


3 4- r 
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I a- as given by typical element (b) of l' 
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APPENDIX C 


AERODYNAMIC EXCITATION PARAMETERS 


C.l GENERAL 

Determination of inter-blade phase angles, excitation reduced 
frequencies, circumferential harmonic indices, and the relation 
between the cosine and sine components of applied airloads are 
described. All cases addressing the turbosystems with flexible 
or rigid hubs/disks placed in uniform or non-uniform inflow are 
considered . 

C.2 INTER-BLADE PHASE ANGLES 

The possible inter-blade angles are determined by factors 
concerning the nature of the imposed aerodynamic excitation . The 
following points are noted here to define the inter-blade phase 
angles: 

1. At any given instant of time, equation (7) of Section 3 
defines the relative velocity components at the leading edge 
point of any chord on the reference/modelled/n=l blade . This 
blade is also identified in Figure 2.5. 

2 . At the same instant , the relative velocity components at the 
corresponding point on any other (say, n ) blade can also 
be obtained from equation (7) , Section 3 , when 
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is replaced by 


where 


a) T 


BL 


-p B L ^ in 

s 


[T^"] = [Tf ] 


0 




o o 

cosy^ 

- siv\ cos yj 


J . 


with 


/ ^ J / / , 


and 


(ci) 


(C7.) 


b) -IwAj-Tunnel and O^A are considered for the n ^ 
blade . 


3. Finally, the ratio of the oscillatory part of velocity 
component VA g- on the nth. blade to that on the (n+1) th. 
blade yields the inter-blade phase angle a~ as 


osc ( va )— h\ocle ^<r' 

: = & 

osc ( v/A)- (n^-O^ bla^e 

^ ) Th. 


(C3) 


This definition of o" is consistent with that of Ref . 12 . 

The above procedure is used in defining cT in each of the 
following cases . 


C.3 RIGID HUB/DISK, UNIFORM INFLOW 


Excitation Reduced Frequency 
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With the turbosystem axis of rotation misaligned at a small 
angle with the uniform inflow velocity, the only excitation 
frequency is given by 

40 = j -O- I (c4-) 

where ~CL is the angular velocity of rotation of the 
turbosystem. 

The corresponding reduced frequency for the blade is defined 
with any of its chords selected as reference: 


k 


bla^la 








(cs-) 


where -^^,,e^and are the semi-chord and cascade relative 

inflow velocity, respectively, for the reference chord . 


Reduced frequencies for other chords are scaled from the above 
blade value . 


Inter-Blade Phase Angle 

This is given by 

<r =: - 27T 

h\acke ' > 


with N being the total number of cyclic segments in the 
structure . 

Circumferential Harmonic Index 
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The only permissible value is 0 , together with the 
requirement that the degrees of freedom on the segment boundaries 
be constrained to zero. 

-kc -ki 

Relation Between r and P 

These components of applied airloads are described in detail 

in Section 2.5. For the rigid hub/disk and uniform inflow case, 

— ks ^kc 

r is identically zero, and only i exists. 


C . 4 RIGID HUB/DISK, NON-UNIFORM INFLOW 


Excitation Reduced Frequency 

The possible excitation frequencies are based on the 
circumferential harmonic contents of the spatially non-uniform 
inflow, and are given by 




where p takes on the positive integer values of the 
above-mentioned inflow harmonics. 


Reduced frequencies are as per equation (C5) with from 

equation (C7) . 

Inter-Blade Phase Angle 

Permissible inter-blade phase angles are obtained from = 
. iTt/fj . However, in order to conveniently determine the 
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associated structural circumferential harmonic index, P 

:: The flexible huh/disX non-uniform inflow case, the inter-hlade 

written as follows: 
phase angles are wrxtu 




b; = c^-o/^ 

_ M /a. 


fsj 

fsi evey\ 


o 


< b ^ b 


a- - 


^ . -xtr/fO 


; . 0 1 , 2-1 
<1 ' / ^ ^ 


(eg) 




t9) 


^ j. j. ^ - __ — — . — — 

• n t-oaether with the requirement 
The only permissible value is 0, together 

of freedom on the segment boundaries be 
that the degrees of freedom 

constrained to zero. 



Only P "" type loads exist 


C.5 vrnVTRLE HUB/DIS Ki-SSI^g^^g^^ 

..citation Reduc ed_Freauei^ 

The only admissible frequency is 

40 -■ I 1 > 
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with the reduced frequency given by equation (C5) . 

Inter-Blade Phase Angle 

Only one inter-blade phase angle exists; 

CT* — — '^TT /fJ 

Circumferential Harmonic Index 

The only possible index can be written as 

k = I cr I / = 1 {cs'L) 

o 

Relation Between r and p 

These two load components are related as 
be 

P = ± i P 

where the + sign applies if o , and the - sign applies 

otherwise . 

C.6 FLEXIBLE HUB/DISK r NON-UNIFORM INFLOW 

Excitation Reduced Frequency 
The possible excitation frequencies are 
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CO - 


I?. I 2i I (CAAr) 

where p is a positive integer reflecting the circumferential 
harmonic contents of the spatially non-uniform inflow. 

Reduced frequencies, in turn, are obtained from equation (C5) . 

Inter-Blade Phase Angle 

The permissible values are computed as per equations (C8) and 
(C9) . 

Circumferential Harmonic Index 

For selected excitation harmonics p, the permissible 
circumferential harmonics representing the structural motion are 
given by 

k. = 1 O' I / ( 2TT /hi) 

/ 

where (T^'s are from equations (C8) and (C9) . 

Tckc — ks 

Relation Between r and P 

This is described by equation (C13) . 
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SYMBOLS 


A 

B 

C 

c 

G 

g 

I 

i 

K 

k 

1 

M 

Ml 

N 

n 

P 

Q 

T 

t 

u 

V 

w 

w 


x,Y,z 


Chordwise aerodynamic matrix for aerodynamic modes 
Damping matrix 

Coriol'is acceleration coefficient matrix 

Matrix of ■modal pressure coef f icients ^ constants 

Chord 

Transformation matrix 
Structural damping 
Chordwise weighting matrix 

Stiffness matrix 

Circumferential harmonic index , reduced frequency 
Semichord 

Mass matrix , cascade relative inflow Mach number 
Centripetal acceleration coefficient matrix 
Number of blades on turbosystem 
n ^ cyclic sector 
Load vector 

Aerodynamic coefficient matrix for structural modes 

Transformation matrix 

Time 

Physical degrees of freedom 
Cascade relative inflow velocity 
Virtual work , weighting matrix 
Strip width, downwash velocity 
Chord local coordinates 
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Y 

1 

A 


9 ^ 

Q 


Inclination angle of turbosystem axis of rotation 
Aerodynamic modal coordinates 
Sweep angle 
Stagger angle 

Structural modal coordinates 
Hass density of flow 
Inter-blade phase angle 
Chordwise aerodynamic modes 
Structural modes 

Pressure distribution in aerodynamic modes 
Rotational speed 
Circular frequency 


9 

i 

P 

s 


Subscripts 

Grid point on chord, chordwise aerodynamic mode 

Structural mode 

Cascade theory pressure point 

Chord s 


B,b 

G,g 

K 

Lrl 

n 

s 


Superscripts 
Basic coordinate system 
Global coordinate system 

Independent solution set in 'symmetric components 
Local coordinate system 

n cyclic sector , normal to chord 

th. 

Chord s 
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Fourier coefficients ( 'symmetric components * ) 


o 


kc 

ks 




N/2 J 
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